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Abstract 

We analyze the beating behavior of the magnetic quantum oscillations in a 
layered metal under the conditions when the cyclotron energy fiw c is compara- 
ble to the interlayer transfer energy t. We find that the positions of the beats 
in the interlayer resistance are considerably shifted from those in the magneti- 
zation oscillations, and predict that the shift is determined by the ratio huj c /t. 
A comparative study of the Shubnikov-de Haas and de Haas-van Alphen ef- 
fects in the quasi-two-dimensional organic metal /3-(BEDT-TTF)2lBr2 ap- 
pears to be consistent with the theoretical prediction. 
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In the last decade, the de Haas- van Alphen (dHvA) and Shubnikov-de Haas (SdH) effects 
were extensively used for studying quasi-two-dimensional (Q2D) organic metals (see for a 
review Due to extremely high anisotropies of the electronic systems the amplitudes 

of the oscillations are strongly enhanced in high-quality samples of these materials (see e.g. 
PH) and often cannot be described by the Lifshitz-Kosevitch (LK) formula derived for 
conventional three-dimensional metals |J. A number of theoretical works performed in the 
last years on the dHvA effect in two-dimensional (2D) and Q2D systems []7|— |T^] provide a 
consistent theory which can be used for a quantitative analysis of experimental data obtained 
on Q2D systems as long as many-body interactions only lead to constant renormalization 
effects and magnetic breakdown effects are not concerned. 

The situation with the Q2D SdH effect is more complicated; here even some qualitative 
questions remain open. One of them is the origin of the phase shift in the beats of the 
resistivity oscillations with respect to those in the magnetization. The beating behavior 
of the oscillations in Q2D metals is well known to originate from the slight warping of 
their Fermi surfaces in the direction normal to the 2D plane. The superposition of the 
contributions from the maximum and minimum cyclotron orbits is expected to lead to an 
amplitude modulation of the fc-th harmonic by the factor cos(2irkAF/2B — 7r/4), where 
B is the magnetic field and AF = (ch/2ne)(A max — A min ) is the difference between the 
oscillation frequencies caused by the extremal orbits with the fc-space areas A max and A min , 
respectively || . From the beat frequency one can readily evaluate the warping of the Fermi 
surface (FS) and hence the interlayer transfer integral, At ~ sfAF/F (see e.g. [3>0])- The 



situation becomes less clear when the warping is so weak that less than one half of the beat 
period can be observed experimentally. In principle, an observation of one single node would 



already be quite informative [14], provided the phase offset (i.e. the phase of the beat at 



1/B — > 0) is known. In the standard LK theory this phase offset is strictly determined by 
geometrical reasons to be equal to — 7r/4 for both dHvA and SdH effects ||. 

However recent experiments on layered organic metals k-(BEDT-TTF)2 Cu[N(CN)2]Br 
14| and (BEDT-TTF)4[Ni(dto)2] |TjJ have revealed a significant difference in the node po- 



sitions of beating dHvA and SdH signals. The respective phase shift in the latter compound 
was estimated to be as big as ir/2. Noteworthy, in both cases the oscillation spectrum was 
strongly dominated by the first harmonic when no substantial deviations from the standard 
LK theory is expected. 



Although several potential reasons for this behavior have been outlined in Refs. [|14}|T5 



the most plausible one seems to be its association with the Q2D nature of the electronic 



system |JL4}| . However, the limited amount of the reported experimental data is not sufficient 
for a detailed analysis. Moreover, multiply connected FSs typical of both compounds might 
lead to additional complications due to effects of magnetic breakdown, interband scattering 
etc. 

In order to clarify the problem, we have carried out comparative studies of the oscillat- 
ing magnetization and interlayer resistivity of the radical cation salt /3-(BEDT-TTF) 2 IBr 2 . 
This material exhibits a relatively simple behavior without superstructure transitions or 
insulating instabilities. It is normal metallic down to low temperatures and undergoes a su- 
perconducting transition at T c pa 2 K [Tj| . Its electronic properties are basically determined 
by a single cylindrical FS |L7j slightly (by ~ 1%) warped in the direction perpendicular 
to the highly conducting BEDT-TTF layers [|T^,^IJ. Thus, the present compound appears 
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to be an ideal object for our purposes. The experiment was done on a high-quality single 
crystal /3-(BEDT-TTF) 2 IBr 2 at T « 0.6 K in magnetic field up to 16 T. To assure exactly 
the same conditions for the dHvA and SdH effects (in particular, identical field orientations 
are of crucial importance for our purposes!) the measurements were performed in a set-up 
providing a simultaneous registration of the magnetic torque and resistance |H| . In the field 
range between 7 and 16 T we have observed clear beating with several nodes in both dHvA 
and SdH signals. The positions of the beat nodes in the SdH signal are found to be different 
from those in the dHvA signal, the difference being dependent on the magnetic field. 

Below we propose an explanation of the phenomenon based on the consideration of 
both density of states (DoS) and Fermi velocity oscillations contributing to the interlayer 
magnetotransport in a Q2D metal and then compare the theoretical estimations with the 
experimental results. 

We consider a Q2D metal in a magnetic field perpendicular to the conducting layers with 
the energy spectrum 

£n,k z — h^c { n + 1/2) — 2t cos(k z d) (1) 

where t is the interlayer transfer integral, k z is the wavevector perpendicular to the layers, 
d is the interlayer distance, uj c = eB/m*c is the cyclotron frequency. Both hu c and t are 
assumed to be much smaller than the Fermi energy. 

The DoS of electron gas with this spectrum can be easily obtained performing the sum- 
mation over all quantum numbers at a fixed energy: 

9(e) = £ = (2) 

«=o yV - (e - huo c (n + 1/2)) 2 

where Nll is the Landau level degeneracy. The sum over Landau levels can be represented 
as a harmonic series using the Poisson summation formula [^]. As a result one gets |11[| : 



,( e) ocl + 2E(-l) fe cos^ 7 -jjo( 7 -j (3) 



k=l 



We shall consider the case 47rt > tkv c when the beats can be observed. Then the zeroth 
order Bessel function J (nk4t/huj c ) describing the beating of the DoS oscillations can be 



simplified, as Jq(x) ~ y 2/7rxcos(x — tt/4). Further, we consider the limit of strong harmonic 
damping, retaining only the zeroth and first harmonics. In this limit the oscillations of the 
chemical potential /i can be neglected. Knowing the DoS we can now obtain the oscillating 



part of the magnetization as [[Tl[ 

cos (P- ~ ~) r t ( 4 ) 
\hujj \hu c a) 1 K ' 



M oc sin 



where Rt is the usual temperature smearing factor. This expression coincides with the 
result of the three-dimensional LK theory 0] and allows an evaluation of t from the beat 
frequency. 

The interlayer conductivity o zz can be approximately evaluated from the Boltzmann 
transport equation |^T|, assuming the impurities are point-like: 
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o zz = e 2 Yj v 2 z (k z )5(e(n,k z ) - fi) • r(/i) = e 2 /( / u)r(/i) 

m={n,k z ,k x } 



(5) 



where v z is the z-component of the electron velocity and r(/z) the momentum relaxation 
time at the Fermi level. The latter, in Born approximation, is inversely proportional to the 
DoS: l/r(/i) oc g(fi) and oscillates in magnetic field according to Eq.(||]). In addition, when 
the cyclotron energy is comparable to the warping of the FS, the oscillations of the electron 



velocity summed over the states at the Fermi level, I(n) = 



become also important 



To calculate this quantity one has to perform the integrations over k x and k z in (|5|) and 
substitute the expression for v z : 



; (e,n) = ^ V /4t 2 -(e-^ c (n + l/2)) i 



(6) 



As a result one obtains 



Tl , ^ N LL d 2 r—— — - 777^1 ±N LL dH 

1(e) = 2^ n _ t y At 2 — (e — huj c (n + 1/2)) = — x 



n=0 



2nh 



4vrf 
8hu r 



+ E 

k=l 



2k 



cos 



'2irke\ ( Ankt" 



TlLOr 



(7) 



Here we again applied the Poisson summation formula |22]. 

The first order Bessel function Ji(4irkt/hLJ c ) entering (0) also describes beatings but with 
a phase different from that of the DoS beatings given by J Q (4nkt/huj c ) in Eq.@: at large 
x, J\(x) ~ ^2/Tixsm(x — 7r/4) is just shifted by tc/2 with respect to Jo(x). 

The phase shift in the beating of the oscillations in g(e) and /(e) is illustrated in Fig. 
1 in which these quantities are plotted for two different values of the ratio At/huj c . When 
At/huj c = 2.25 the DoS oscillations have a maximum amplitude of the first harmonic (Fig. 
la). At the same time the oscillations of 1(e) exhibit a nearly zero amplitude of the first 
harmonic as shown on Fig. lc whilst their second harmonic (not shown in the figure) is at 
the maximum. By contrast, when 4t/hu c = 1.8 (Fig. lb,d) the first harmonic of the DoS is 
at the node whereas that of 1(e) has the maximum amplitude. 

The difference between the phases of the beats in oscillating g(e) and /(e) leads to a 
shift of the beat phase of the SdH oscillations with respect to that of the dHvA oscillations. 
Indeed, taking into account that r oc 1/g, substituting Eqs. (^) and (|7|) into (|^), applying 
the large argument expansions of Jo(x) and J\(x) and introducing the temperature smearing, 
we come to the following expression for the first harmonic of the interlayer conductivity: 



o 7 



/27T/i\ /47Tt 

oc cos - — COS 



7T 

4 + 



Ri 



(8) 



where 



arctan(a) and a = hu c /2irt 



(9) 
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Comparing these expressions with Eq.(f|), one can see that the beats in the SdH and dHvA 
oscillations can become considerably shifted with respect to each other as the cyclotron 
frequency approaches the value of the interlayer transfer integral. 

The above evaluation is based on the semi-classical Boltzmann equation and certainly 
is not expected to give a precise result. Nevertheless, as will be seen below, its predictions 
concerning the phase shift are in good qualitative agreement with the experiment and we 
therefore believe that it correctly reflects the physics of the phenomenon. 

Fig. 2 shows the oscillating parts of the magnetization and interlayer magnetoresistance 
in the normal state of /5-(BEDT-TTF) 2 IBr 2 at magnetic field tilted by 9 ~ 14.8° from the 
normal to the BEDT-TTF layers. The curves have been obtained by subtracting slowly 
varying backgrounds from the measured magnetic torque t(B) and resistance R{B) and 
(for the magnetization oscillations, M oc r/B) subsequent dividing by B. The fast Fourier 
transformation (FFT) spectra shown in the insets reveal the fundamental frequency of ~ 
3930 T in agreement with previous works [|T^,^|. The second harmonic contribution is 
about 1% of that from the fundamental one at the highest field. 

Clear beats with four nodes (indicated by arrows) are seen in both the dHvA and SdH 
curves. We have assured that the observed beats originate from the warping of the cylindrical 



FS by checking the angular dependence of their frequency |2"0"|. The fact that the oscillation 



amplitude does not exactly vanish at the nodes was attributed to slightly different cyclotron 



masses at the extremal orbits of the FS [20|. The positions of the nodes determined as 



midpoints of narrow field intervals at which the oscillations inverse the phase are plotted on 
Fig. 3. The straight line is a linear fit of the magnetization data revealing the beat frequency 
AF = 40.9 T that, according to Eq.(^), corresponds to At/eF = AF/F = 1/96. The error 
bars in the node positions do not exceed ±3 x 10~ 4 T _1 for N = 3 to 5 and are somewhat 
bigger, ~ ±10~ 3 T _1 , for N = 6 due to a lower signal-to-noise ratio. We note that although 
the angle 9 = 14.8° corresponds to a region in the vicinity of the maximum beat frequency 
(ma.x{AF(9)} « 42.0 T for the given field rotation plane), the sensitivity of the node 
positions to the field orientation is still quite high: the nodes shift by ~ 2.3 x 10 -3 T _1 at 
changing 9 by 1°. Thus, if one has to remount the sample between the torque and resistance 
measurements, even a slight misalignment may cause a substantial additional error. In our 
experiment both quantities were measured at the same field sweep, hence, such an error was 
eliminated. 

From Figs. 2 and 3 one can see that the nodes of the SdH oscillations are considerably 
shifted to higher fields with respect to those of the dHvA oscillations. The shift grows with 
increasing the field. Both these observations are fully consistent with the above theoretical 
prediction. 

In order to make a further comparison between the experiment and theory, we plot the 
quantity tan(</>) (where <ft is the phase shift between the beating of the SdH and dHvA 
oscillations obtained from Fig. 3) as a function of magnetic field in Fig. 4. A linear fit to 
this plot (dashed line in Fig. 4) has a slope of 0.037 1/T |p3| . A substitution of this value 



and the cyclotron mass m* = 4.2m e , obtained from the temperature dependent amplitude 
of the fundamental harmonic, into Eq.@ yields an estimation for the interlayer bandwidth 
At ~ 0.48 meV or the ratio At/ep = AF/F ~ 1/230. This is somewhat smaller than the 
value 1/96 obtained directly from the ratio between the beating and fundamental frequencies. 
However, taking into account an approximate character of the presented theoretical model, 
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the difference is not surprising. Further theoretical work is needed in order to provide a 
more explicit basis for the quantitative description of the phenomenon. 

Summarizing, the beats of the SdH oscillations in /3-(BEDT-TTF) 2 IBr 2 are found to be 
shifted towards higher fields with respect to those of the dHvA signal. We attribute this 
effect to interfering contributions from oscillating DoS and Fermi velocity to the interlayer 
conductivity of this layered compound. Thus, the observed behavior appears to be a general 
feature of Q2D metals which should be taken into account whenever the cyclotron energy 
becomes comparable to the interlayer transfer energy. 

We are thankful to A.M. Dyugaev, I. Vagner and A.E. Kovalev for stimulating discus- 
sions. The work was supported by the EU ICN contract HPRI-CT- 1999-400 13, and grants 
DFG-RFBR No. 436 RUS 113/592 and RFBR No. 00-02-17729a. 
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Figure captions 



Fig. 1. a - the DoS (solid line) near the Fermi level and its first harmonic (dashed line) 
according to Eq.([|), at the ratio 4t/hu c = 2.25; b - the same at the ratio 4t/hu c = 1.8; c 
- the quantity /(e) (solid line) and its first harmonic (dashed line) according to Eq.(7), at 
At/huj c = 2.25; c - the same at 4t/hu c = 1.8. In all four panels: the dotted lines are the 
contributions from individual Landau levels. 

Fig. 2. dHvA (left scale) and SdH (right scale) oscillations in /3-(BEDT-TTF) 2 IBr 2 at 
9 14.8°. Insets: corresponding FFT spectra. 

Fig. 3. The positions of the nodes in the oscillating magnetization (filled symbols) 
and resistance (open symbols) versus inverse field. The straight line is the linear fit to the 
magnetization data. 

Fig. 4. Tangent of the phase shift between the node positions in the SdH and dHvA 
signals taken from the data on Fig. 3 as a function of magnetic field. The dashed line is a 
linear fit according to Eq.(|5|) P5| . 
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(a) At = 2.25Jkb c 
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Fig. 1 of ..Anomalous beating phase..." by P. Grigoriev et al. 




Fig. 2 of "Anomalous beating phase..." by P. Grigoriev et al. 




Fig. 3 of "Anomalous beating phase..." by P. Grigoriev et al. 
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Fig. 4 of "Anomalous beating phase..." by P. Grigoriev et al. 



